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Abstract 

We present a study of diffusion of small tagged particles in a solvent, 
using mode coupling theory (MCT) analysis and computer simulations. The 
study is carried out for various interaction potentials. For the first time, us- 
ing MCT, it is shown that for strongly attractive interaction potential with 
soft core (allowing interpenetration between the solute-solvent pair) the dif- 
fusion exhibits a non-monotonic size dependence. This was earlier predicted 
in simulation and experimental studies and was connected to levitation effect 
[J. Phys. Chem. B 2005, 109, 5824-5835]. Our MCT analysis reveals that for 
weak or no attractive interactions, all the small solute particles studied here 
show levitation through the inter-solvent transient cage. However, for strong 
attractive interaction the levitation is not present for the smallest particle 
sizes. It is found that for systems where the interaction potential is hard, 
not allowing any interpenetration between the solute-solvent pair, the solute 
cannot explore the inter-solvent cage. Thus these systems will never show 
any non monotonic size dependence of diffusion. We also show that although 
levitation is a dynamic phenomena, the effect of levitation can be obtained 
in the radial distribution function. 
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1 Introduction 



Diffusion phenomena in ffuid has immense importance in various fields such 
as industrial chemistry, soft matter physics, and biology. The motion of 
tagged particles in various mediums have been extensively studied using, mi- 
croscopic theories, simulations and experiments [1-12]. These studies have 
predicted many interesting results. The theoretical studies of diffusion using 
mode coupling theory have predicted that the microscopic origin of many 
phenomena are quite different from what it had been envisaged earlier. For 
example, the diffusion in a neat liquid (where the solute particle is one of 
the solvent molecules) is usually found to follow the Stoke-Einstein relation 
(SER) (Docl/R where R is the radius of the particle). However, the SER 
should ideally be valid in the hydrodynamic regime where the solute particle 
is so big that it perceives the solvent as a continuum without any microscopic 
details. It was found that SER is valid because the viscosity and the friction 
are determined by the same dynamical variable, namely the density auto 
correlation function [I] [2]. Thus the study revealed that the validity of SER 
has a microscopic rather than a macroscopic hydrodynamic origin. The dif- 
fusion of small solute particles in a sea of larger solvents are known to violate 
the SER. Experiments have repeatedly shown that for small solute particles 
the SER significantly underestimates the diffusion coefficient [21 S] and this 
has been attributed to the micro viscosity effect or fractional viscosity ef- 
fect. Analytical MCT studies have shown that the decoupling of the solute 
motion from the solvent dynamics is the origin of the enhanced diffusion 
[5]. The diffusion of small solutes not only show the breakdown of the SER, 
but also other anomalies. There has been a series of studies of diffusion of 
small solutes in different mediums by Yashonath and co-workers [B] [Z] [S] . The 
medium in which diffusion was studied varied over a wide range from zeolites 
to other porous solids and also included Lennard- Jones (LJ) solvent. The 
simulation studies revealed that diffusion of small solutes show a diffusivity 
maxima, independent of the nature of the medium. This counter intuitive 
behaviour of solute diffusion through different types of porous structures was 
first reported by Henson et al. [S]. Yashonath et al. have systematically 
varied the size of the solute by keeping the porous zeolite structure the same 
[6]. In their study of diffusion through zeolite Y they have found that for a 
spherical sorbet of diameter 6.0^4° which has dimensions comparable to the 
bottleneck of the zeolite Y, the force exerted by the zeolite on the diffusing 
particle is minimum which leads to the diffusion maxima. In general they 
noticed that in zeolites when the solute size is about 80% of the size of the 
neck of the zeolite the diffusion is maximum. In their study of solute diffu- 
sion in LJ solvent they found that when the solute size is about 0.24 times 
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that of the solvent the diffusion shows a maxima [Qj. From Voronoi Poly- 
hedra analysis they have obtained a void and neck size distribution of the 
inter-solvent transient cages and have shown that, when the solute-solvent 
diameter is smaller but similar to the size of the neck, the solute can dif- 
fuse faster. They have attributed this diffusivity maxima to levitation effect 
(LE), stating that when the solute size and the size of the neck are similar 
then the attractive force felt by the solute while passing through the pore is 
uniform from all directions and this helps the particle levitate through the 
pore. Whereas, in case of smaller solutes, proximity to one part of the neck 
leads to a large asymmetry in the force felt by the solute. The solute tends 
to spend longer time bound to one side of the neck or one solvent particle (in 
case of LJ solvent), which leads to lower diffusion constant. Lynden Bell and 
co-workers have also found such diffusivity maxima for ion diffusion in polar 
solvent (water) [11]. Experiments have also found this type of non-monotonic 
size dependence of diffusion [9]. All these studies show that the non mono- 
tonic size dependence of diffusion is a generic phenomena independent of the 
nature of the solvent. 

To understand this behaviour observed in simulation and experimental 
studies we require a detailed knowledge of the diffusion mechanism. Earlier 
studies have shown that mode coupling theory (MCT) can successfully ex- 
plain diffusive dynamics and its anomalies over a wide range of solute-solvent 
size ratios pQ. Hence it is imperative to study this non monotonic size depen- 
dence of diffusion and its connection to levitation, using the mode coupling 
theory. 

In this article we use the fully self-consistent microscopic mode-coupling 
theory and molecular dynamic simulations to understand the origin of non- 
monotonic size dependence of diffusion and its relationship with levitation. 
Although the present study is confined only to solute diffusion through pure 
solvents, we expect that the understanding of the diffusion mechanism in this 
system will help us in future to study the solute diffusion in porous medium. 
The present analysis involves four different systems by systematically varying 
the solute diameter for each system. The interactions between the solute - 
solvent and solvent-solvent pairs are modeled in terms of LJ potential, where 
the interaction parameters are varied from large attractive to repulsive values. 
The solute-solvent radius ai2 is varied in such a way that in some cases inter- 
penetration between the solute-solvent pair is allowed. We show that for 
repulsive and weak attractive interaction the LE is present for all the solute 
sizes, whereas for large attraction only certain solutes of intermediate sizes 
can levitate through the transient inter-solvent cage. The smaller solute 
particles cannot levitate and have larger probability to primarily stick to 
one solvent molecule. The study also reveals that although levitation is a 
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dynamic phenomena, its signature is present in the static radial distribution 
function (RDF). 

The work is arranged as follows. The simulation details are given in 
Sec. 2. Sec. 3 contains theoretical analysis. In Sec. 4 the results are presented 
and finally, Sec. 5 concludes with short discussion about the result. 



2 Simulation Details 

The atomistic models which are simulated are two component mixtures where 
particles of type i interact with those of type j with pair potentials, {Uij{r)) 
where r is the distance between the pairs. The potentials are shifted and 
truncated Lennard- Jones potentials. 



where u[^'\r- a^,, e,,) = 4e,, [(2^^)12 _ (fii)6] = 2^ a,, for Week Chandler 



Anderson (WCA) system [22] and rf^ = 3.5aij for LJ system. Were i,j=l 



for solvent and i,j=2 for solute. 

Four systems, consisting 456 solvent or host atoms {Nh) and 44 solute or 
guest atoms (Ng) are studied. + Ng = 500. 

The different models are distinguished by different choices of lengths and 
energy parameters. Lengths, temperature and time are given in units of an, 
and respectively. The details of the models are given in the 

Table (I) and mi = m2 = m = 1 for all the cases. Note that for system 
1, 2 and 3, the interaction between the solute and the solvent has a soft 
core which allows inter-penetration between solute solvent pair. However in 
system 4 the solute-solvent diameter is given by Lorentz-Berthelot (LB) rule 
[23] 

If eii=0.25 KJ/mole,and (Tii=4.1 A^, (T22 is varied between 0.3 to 1.5 
A' ,mi = 40amM, then r=5.186 ps and the system is identical with set-4 
studied by Yashonath et.al [B]. 

The sample is cubic with side length 7.86789(Tii. We have carried out MD 
simulations at reduced temperature, T*=1.663, reduced density for the sol- 
vent, p*=0.933,in micro-canonical ensemble(NVE) using velocity integration 
with time step of O.OOlr. 

For all the simulations mentioned above,the system has been equilibrated 
for 1 ns followed by production run of 2 ns. The molecular dynamics simu- 
lations carried out using LAMMPS package [1] 
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Tabic 1: Parameters for model systems used in this study, where (Th, (T22 and (T12 are the 
solvent, solute and solute-solvent diameter respectively and en, £22 and ei2 are the interac- 
tion parameters for the solvent-solvent, solute-solute and solute-solvent pairs, respectively. 
For system 1 , 2 and 3 inter-penetration between the solute-solvent pair is allowed and sys- 
tem 4 follows LB rule for mixing and the 0^2 values are 0.073, 0.098, 0.122, 0.171, 0.200, 
0.220, 0.244, 0.293, 0.317, 0.366. 
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The simulation results have also been checked for Ng=l and the results 
(like RDF and D) remain identical with that for A''g=44. 

3 Theoretical Analysis 

For the mode-coupling calculation we consider a single solute (guest) particle 
in a sea of solvent (host) particles. The diffusion coefficient (D) of the tagged 
solute is given by the well known Einstein relation 

D = kBT/mV{z = 0) (2) 

Where T{z) is the frequency dependent friction on the particle having mass 
'm'. The frequency dependent friction T{z) can be calculated using mode 
coupling theory. The separation of time scale between the binary and re- 
peated recoUisions are invoked to decompose the friction into short and long 
time part. The short time part arises due to direct binary collisions, and a 
long time part which arises due to correlated recollisions. The calculation of 
the recoUision term is non trivial and includes information of the five hydro- 
dynamic modes. It has been reported earlier that for small solute particles, 
the density mode of the solvent primarily contributes to the recollisional fric- 
tion. If we neglect the contribution from other hydrodynamic modes then 
the friction can be written as, [S] 

r^{z)c^T%{z)+Rj^{z) (3) 

where r^(z) is binary collisional friction. R'^p{z) arises from the coupling of 
the tagged particle motion with the density fluctuation of the medium. 

A single binary collision between the tagged particle and a solvent particle 
, in presence of the other solvent particles is described by F^. The calculation 
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of the full time dependence of is non trivial. However it has been shown 
that the re collision term begins as [T3], thus the binary collision term 
contains all the contributions to the order t^. As only even powers of t 
appears thus the binary collision can be assumed to be a Gaussian, 

m) = <.exp(-t7r./) (4) 
here cjoij is the Einstein frequency : 

= 0) = u^l^ = £^l dvg,,{r)V\,{T) (5) 

where %(r) is the inter-particle potential and gij{r) is the static pair correla- 
tion function (RDF) between the i and the j type particles, p is the particle 
number density. In Eq. [5] the relaxation time Tij is determined by taking 
the second derivative of a more general expression for binary friction and is 
given by [T5] . 

+ (l/6p) J [dq/(2vr)3]7g(q)(5ii(q) " ^h^M) (6) 

where summation over repeated indices is implied. 5'ii(q) is the static struc- 
ture factor for the pure solvent and ■y'^f- is given by 



^diM) = -{p/^) j rfrexp(-zq.r)^ij(r)V"V''%(r) 

= rfl'M) + - rfWM) (7) 
where 7^i^.(q) = Ydtji^) and 7^i^.(q) = 7;^^.(q) 

The term arising in the friction due to the coupling to the density fluctu- 
ations in the surrounding medium is given by 

R%it) = (p/m/3) J [rfq7(2vr)3](q.q')'g'^ii(g) 

X h,{q')]'[<j>'{q',mq\t) - t)0o(g', t)] (8) 

Where (3=-^^ and q' is an unit vector along the arbitrary z direction and 
Cij{q) is the two particle direct correlation function between particles i and 
j. (f){q,t) and (j)'^{q^t) are the full and the self intermediate scattering func- 
tions respectively. 0o(g, t)=exp(- ^^ff^^^^ ), and (pl{q,t) = exp{-^) are 
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the inertial part of the intermediate scattering function and it's self part, 
respectively. 

(psiQyt) and 4>{q,t) are calculated self consistently using the well known 
MCT equations (^^. 

0,(t) + r^n(g,2 = o)0,(t) + fi2^,(t)+ / M,,j{t - T)^g{T)dT = (9) 

Jo 

0%(t)+r\(g,^ = O)0^,(t)+fi^0%(t)+ f Mlijit - T)<P%{T)dr = (10) 

Jo 

where = ji^k^T ^2 _ g^MZ and Ma^i(t - r) and M'.At - r) are 
the memory kernels. They give rise to the long-time tails in the intermediate 
scattering function and its self part. The short time part is considered to 
be delta-correlated with the strength given by Tfj{q,z = 0). r^(q,z=0) is 
written as 

rg(g, ^ = 0) = / rft<exp(-tVr5(g)) (n) 

Where Tjj(g) is calculated from 

+ (l/6p) J [rfq/(2vr)3]7g(q)(5ii(q)-l)7g(q) (12) 

MMt) = ^|rfq[k-(qcn(g) + (k-q)cn(fc-g))f 5n(g)5n(fc-g) 
[0(g, t)<P{q -k,t)- Mq, t)M(l - k, t)] (13) 

Ml,,{t) = ^ I dqSn{q)cM\^?[M<l - k,m{q,t) - 0^(g - k,t)M<l,t)] 

^ (14) 
Eq (fT3|) is the memory function for the {/)(g, t) calculation and Eq fll4l) is 
for (j)s{q,t) 

The integration in the memory-functions are converted to Riemann sums 
for the discrete momentum values [13]. We take 500 grid points with q 
starting from 0.2 upto a cutoff of 199.8 with step size dq=0.4. Usually MCT 
calculations are done with q cutoff of 39.8 but for our calculation since the 
solute sizes are small we need to go until q = 199.8. 

We use gii{r) and gi2{r) from theory and simulation and 5*11 (g) and Cii(g) 
and Ci2(g) are calculated analytically using HMSA closer [16]. We get 4>{q,t) 
and (j)s{<l,t) by solving Eqs.fl9]and [T0|) by self consistent method. 
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4 Result and Discussion 



The Mode coupling Theory (MCT) formahsm discussed above has been 
employed to calculate the self diffusion constants, D. The diffusion value 
for a pure solvent using this fully self consistent MCT at T* = 1.663 and 
p* = 0.933, has been found to be Dmct = 0.194 x 10~^cm^/s which is close 
to the simulated value Dgim = 0.195 x 10~'*cm^/s. 
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Figure 1: Scaled diffusion coefficients obtained from fully self consistent MCT, Dmct 
(circles) and from simulations, -Dsim (triangles) are plotted against (T22 for system 1. Both 
show non-monotonic size dependence of diffusion with a diffusion maxima. To plot both in 
the same graph, the diffusion values are scaled by arbitrary constants. The CT22 values used 
for the calculation are 0.073,0.098,0.122,0.171,0.200,0.220,0.244, 0.293,0.317 and 0.366 



In FigJUwe present both simulated and calculated (using full self-consistent 
MCT) diffusion values against the solute diameter for system 1 (see Table I). 
The studies show a non monotonic size dependence of diffusion. The theory 
and simulation studies both predict that as the solute size is decreased the 
diffusion value first increases and then decreases. This diffusivity maxima as 
discussed earlier has been attributed to levitation effect [6]. Thus it is found 
that MCT can predict levitation effect and diffusivity maxima as found in 
simulations [0] and experimental [12] studies. 

Although the nature of diffusion as a function of the solute diameter, pre- 
dicted by the fully self consistent MCT, is similar to that found in simulation 
studies, the value of diffusion obtained from this method is found to be much 
smaller than that obtained from simulations. Note that in FigJTJthe Dgi^ and 
Dmct are scaled by arbitrary values. Thus, although the diffusion value for 
a pure solvent obtained using fully self consistent MCT, matches with that 
of the simulated results, for smaller solutes the fully self-consistent MCT 
seems to overestimate the friction. Compared to the simulation values, the 
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self intermediate scattering function, (f)s{q,t) obtained using this fully self- 
consistent method (from EqslH] and [TU]) is found to have a longtime tail. For 
smaller solutes this discrepancy between simulated and analytical 0* increase 
in the long time. This implies that the wavenumber dependent friction, -M^ij, 
which is the memory function for 0* and determines its longtime tail is not 
the correct approximation in the limit of small solute particles. There should 
be other terms in the memory function which should lower the value of the 
wavenumber dependent friction, -Mqij, and thus the long time tail. 

To avoid this problem we have calculated approximate 0* using simulated 
mean square displacement < Ar^(t) > and non Gaussian parameter a2 = 



3<Ar4(t)> 
5<Ar2(t)>2 



— 1. The expression is given by [2 



(f)'{t) = exp 



g2 < Ar2(t) > 
6 



^ 1 < Ar^(t) > 

1 + -a2 — — 

2 ' V 6 



(15) 
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Figure 2: Diffusion coefficients obtained from MCT where (j>s{q,t) is calculated using 
Eq |151 DYjc'y' (diamonds), and from simulations, Dgim (triangles) are plotted against 
(722 for system 1. We also plot the bare diffusion coefficient Db (stars). Diffusion values 
are scaled by 10"'* cm^/s. The solute diameters used for this calculation are same as in 
FiglH The diffusion values show non monotonic size dependence. 



With this approximation for the self intermediate scattering function the 
diffusion/friction has been calculated using Eqs. 2-8. Note that the dynamic 
structure factor is always calculated self-consistently. In the calculation the 
simulated radial distribution function has been used as the effect of levitation 
is stronger in the simulated RDF. This point will be discussed later. In FigjS] 
the diffusion values calculated using approximate </)^(t)(from Eq JTSj) . D^/cr^' 
are plotted against the solute diameter for system 1. For comparison the 
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simulated values, Dgi^ are also plotted. Dapprox values are similar to Dgim 
values and does show a non monotonic size dependence. However, the peak 
value is little shifted and for the two smallest sizes (0.073 and 0.098) the 
diffusion values are smaller than that predicted in simulation studies. In the 
same figure the bare diffusion coefficients obtained from inverting the binary 
frictions, {Db = are also plotted. It is found that for small solute sizes 
the binary contribution to friction is dominant and the contribution from 
the density term i?*^^ is negligible. The frictional contributions from the 
other hydrodynamic modes like the longitudinal and transverse current [19] 
have also been calculated and it was found that non of the hydrodynamic 
modes actually contribute to the friction/ diffusion. For small solutes as 
found earlier [5], the solute motion is decoupled from the solvent dynamics 
and the hydrodynamic modes of the solvent do not play any role in the 
solute dynamics. For these small sizes there may be some contribution from 
non-hydrodynamic modes which assist in the particle diffusion. 

Except for the smallest two sizes the behaviour of diffusion constant as 
predicted by the analytical MCT calculation is similar to that predicted by 
the simulation studies. As mentioned in the Introduction, this anomalous 
size dependence of diffusion has been attributed to levitation effect [6] where 
the solute particles with intermediate sizes are expected to levitate through 
the inter-solvent transient pores (cages) whereas the smaller particles are 
expected to be strongly attracted to one solvent particle. 
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Figure 3: Radial distribution functions from simulation are plotted for different solute 
sizes where the solute sizes are varied as given in FiglTJ In the inset the same is plotted 
as obtained using HMSA closer [16]. For both the sets of RDFs the peak value of the 
maxima show a non-monotonic size dependence. The effect is stronger in the simulated 
RDF. 



The effect of the above mentioned phenomena is reflected in the solute- 
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solvent radial distribution function. In FigJHlit is shown that the peak value 
of the RDF does show a non monotonic size dependence of diffusion as a 
function of solute size. The RDFs for the intermediate solute sizes have a 
lower and broader diffused 1st peak, whereas the value of the 1st peak of the 
RDFs for smaller solutes are higher and the peaks are also much sharper. 
We also note the the non monotonicity is more prominent in the simulated 
RDF. However, the analytical RDF also shows a weak non monotonic be- 
haviour thus predicting that the HMSA closer can show a weak effect of the 
levitation. The effect of levitation on RDF will be discussed in details later 
in this section. As discussed earlier, for this range of solute sizes, the bi- 
nary contribution to the friction is dominant. The binary friction (Eq. 4-7) 
is primarily determined by the RDF. Thus the effect of levitation which is 
manifested in the RDF does finally give rise to the diffusion anomaly. 

For system 1 where the anomalous non-monotonic size dependence of 
diffusion is observed, the particles interact via LJ potential and there is a 
strong attraction between the solute-solvent pair (ei2 = 6.0). Note that the 
solute-solvent diameters in this system are smaller than that predicted by 
Lorentz-Berthelot rule [23] • Thus for this system the interaction potentials 
have soft core and interpenetration between the solute and the solvent pairs 
are allowed. 

To find out the individual roles that core softening and strong attraction 
play in the anomalous behaviour of the diffusion, we further study three more 
systems(see Table I). For system 1, 2 and 3 the core softening is allowed, 
however the attractive part of the interaction for system 2 is weaker than 
that for system 1 and for system 3 the interaction is kept only repulsive. 
Between system 1 and 4 , both have strong attractive interaction but system 
4 does not allow any core softening. 

The diffusion values, Dmct, D'^mct' ^ Dgim and Db for systems 2 and 3 
are shown in FigJH Neither simulation nor theory show any non monotonic 
size dependence of diffusion. In both these two systems the smaller particles 
diffuse faster. Just like the diffusion coefficients, the 1st peak values of the 
radial distribution function in Figj5] also does not show any non monotonic 
behaviour. For system 2 and 3, similar to system 1 the friction is primarily 
determined by the binary collisions, which in turn depends on the RDF. Thus 
monotonic RDF implies monotonic diffusion values. From this analysis it can 
be concluded that only in presence of strong attractive interaction between 
the solute-solvent pair, the diffusion shows a maxima for intermediate sizes. 

In order to understand the role of core softening in this anomalous be- 
haviour of the diffusion we simulate a system where the diameter for the 
solute-solvent pair is ai2 = """+°"^^ ^ following Lorentz-Berthelot rule. The in- 
teraction parameter for this system has been kept same as that for system 1, 
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Figure 4: (a) Diffusion coefficients obtained from fully self consistent MCT, Dmct (cir- 
cles), from MCT where (j)s{q,t) is calculated using Eq. [131 ^mct^ (diamonds), and from 
simulations, Dgim (triangles) are plotted against (T22 for system 2. We also plot the bare 
diffusion coefficient Db (stars), (b) Same as in (a) but for system 3. Diffusion values 
are scaled by 10~^cm^ /s. Non of the diffusion values show any non monotonic solute size 
dependence. The solute values are same as in FigHJ 



i.e the solute-solvent interaction has been kept strongly attractive {e^ = 6.0). 
As shown in FiglHti', even though we have strong attraction between the so- 
lute and the solvent, without the core softening the diffusivity maxima is not 
present. The RDF for this system as shown in FigJUb does not show any non 
monotonic behaviour. Unlike in systems 1,2 and 3, in system 4 the friction 
is determined by both the binary and the density contribution terms. 

For the solvent studied here we know that the neck size distribution of the 
transient solvent cage as obtained from the Voronoi Polyhedra analysis has a 
peak around aneck = 0.425crii (see figure 5(b) of Ref. [6]). Yashonath et al. in 
all their studies have defined a levitation parameter 7 = "^tH^ii^ Usually for 
zeolitic system they have found that when the levitation parameter, 7 is close 
to unity the diffusion reaches a maximum value. For solute-solvent system 
they have found that the diffusion value is maximum when 7 0.62. We 
find that for systems which follow LB rule even for the smallest solute sizes, 
the levitation parameter is always greater than unity. Thus these systems 
will never predict anomalous non monotonic size dependence of diffusion. 

In all the systems it is found that Dmct underestimates the diffusion 
value. However, D^j^^^ shows a better agreement with Dsim- For system 
1 this discrepancy is so large that they are not plotted in the same graph. 
We find that the discrepancy reduces with the reduction of the strength 
of attractive force. We also find that the discrepancy is larger for smaller 
solutes. Thus the study reveals that for strong attractive interactions and 
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Figure 5: Radial distribution function (RDF) for system 2 (a) and system 3 (b) for all 
solute particles. The peak values of the RDFs monotonically decrease with the solute size. 
For system 3 which is purely repulsive, the RDF is almost structure less for small sizes. 

smaller solutes the memory function for the self dynamic structure factor 
M^ij as given by Eq. [Tl]is an approximation and should have contribution 
from other terms. 

4.1 Levitation and radial distribution function 

Our analysis reveals that non monotonic size dependence of diffusion is di- 
rectly connected to the non monotonic behaviour of the 1st peak of the RDF. 
To get a deeper insight and also to understand the connection between levita- 
tion and structure of the liquid, we compare the radial distribution functions 
for all the systems. This comparison is done for two solute sizes, (T22 = 0.2439, 
which falls in the normal region where smaller particles always diffuse faster 
and (722 = 0.073 which falls in the region where the smaller particles can have 
lower diffusion values compared to that for larger particles. The idea of this 
analysis is to understand how much the structural information can tell us 
about the dynamics in the system. In the previous subsection it was found 
that for system 2 (where €12 = 1.0) and for system 3 (where the interaction 
between all the particles are purely repulsive) the solute size dependence of 
the diffusion is monotonic. This we believe happens because for these systems 
all the small particles can levitate through the inter solvent cage. However, 
for system 1 only the intermediate size particles can levitate and the smaller 
particles are stuck to one solvent particle. The signature of this can be found 
in FiglTl 

For solute size (T22 = 0.2439, the RDF has a diffused 1st peak. As the 
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Figure 6: (a) Diffusion coefficients obtained from fully self consistent MCT, Dmct (cir- 
cles), from MCT where (l)s{q,t) is calculated using EqUSl -D^f^™^ (diamonds), and from 
simulations, ZJ^im (triangles) are plotted against (T22 for system 4. We also plot the bare 
diffusion coefficient Db (stars). Diffusion values are scaled by 10^*cm^ / s. The diffusion 
values does not show any non monotonic behaviour as a function of solute size. The solute 
diameters arc same as in FiglTJ (b)Radial distribution function (RDF) for system 4 for 
all solute particles. This plot shows well defined RDF with no non monotonic solute size 
dependence of the peak value. 



attraction between the solute and the solvent is increased, for system 1 we do 
see a well defined maxima but a weak signature of the minima is present and 
there is a fiattening of the RDF after the minima. However, for solute size 
0.073, for system 2 and 3 the RDF is almost structure less but for system 1 it 
has a well defined structure with proper maxima and minima. When a solute 
levitates through the inter solvent pore we expect the solute-solvent RDF to 
be structure less which is clearly observed in the RDFs obtained in the WCA 
system (FigJ3)). However, when the solute remains closer to one solvent 
molecule the RDF has a well defined maxima and it cannot levitate through 
the inter-solvent transients cages due to asymmetry in the interactions. The 
value of the maxima directly predicts the percentage of the population which 
cannot levitate. In system 1 although for (J22 = 0.2439 there is a well defined 
maxima but it is weak compared to that for = 0.073. For system 4 where 
LB rule is followed we find that the peak values are almost independent of 
the solute size. If RDF for system 1 and 4 are compared then it is observed 
that for the smaller solute the system 1 has a larger value of the 1st peak and 
for the intermediate solute the system 4 has a larger value of the 1st peak. 
This implies that in system 1, the smaller solute has a larger affinity towards 
the solvent where as the intermediate size solute has a lower affinity. Thus in 
system 1 most of the small particles feel the asymmetric force whereas fewer 
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number of intermediate size particles feel the asymmetric force and a larger 
population of them levitate. 



bO 6 




Figure 7: Solute-solvent radial distribution function (RDF) for four different systems (1,2 
3 and 4), where cr22 = 0.0 73(a) and (722=0.244 (b). For (722 =0.073 the RDF in system 2 
and 3 are diffused but in system 1 and 4 it has well-defined peaks with clear maxima and 
minima. For (722 = 0.244 in system 2 and 3 the peak values are diffused and in comparison 
to the small solute, the peak value for system 1 is much less. For system 4 since there is 
no Icvitation the peak values for both the sizes arc similar. 



5 Conclusion 

In this article we have studied using both analytical MCT and computer 
simulations the size dependence of diffusion of small solute molecules. This 
study is done in order to understand the non monotonic size dependence of 
diffusion of small solutes which was earlier found in simulation studies [S] and 
experiments In the earlier simulation and experimental studies the 

diffusion of small solute particles in porous zeolite structures [I2][26], pure 
solvents[6] and ionic medium[TT] have predicted that the non monotonic size 
dependence of diffusion is a generic phenomena, independent of the nature 
of the solvent. Yashonath and co-workers [6] have connected this behaviour 
to levitation effect, stating that when a solute size is approximately similar 
to the size of the zeolitic pore or in case of pure solvent, the size of the inter 
solvent transient cage then the solute particle feels symmetric attraction 
from all directions and thus levitates through the pore and diffuses faster. 
However, a smaller solute gets closer to one wall of the pore or one solvent 
molecule and thus feels an asymmetric force and spends a longer time in the 



15 



cage before diffusing out. This leads to the lower diffusion value for a smaller 
solute. 

In order to understand this phenomena we use both analytical MCT and 
simulation techniques. Study has been done for 4 different systems. In three 
systems (1, 2 and 3) the solute-solvent diameter is such chosen that interpen- 
etration between the solute and the solvent is allowed. It is found that for 
system 1 which has strong attractive interaction the diffusion coefficient does 
show a non-monotonic size dependence of diffusion. For system 2, which has 
weak attractive interaction and system 3 which has only repulsive interaction 
there is no non-monotonic behaviour as obtained from both analytical and 
simulation studies. This implies that systems with only strong attractive 
interactions show this anomalous behaviour. 

To understand the effect of core softening, we have simulated a system 
which is strongly attractive but the solute-solvent diameter is calculated us- 
ing LB rule. In this system again we do not find any non-monotonic size 
dependence. Note that, for a solute to levitate, the levitation parameter 7 
needs to be close but less than unity [25]. For any system where LB rule 
is followed the 7 value can never be less than unity as the solute-solvent 
diameters are always bigger than the inter-solvent cage. Thus we conclude 
that for solute-solvent system to observe this anomalous size dependence the 
attractive interaction between the solute and the solvent should be strong 
and also interpenetration between the solute-solvent pair should be allowed 
for the solute to explore the inter-solvent cage. This core softened potential 
or interpenetration should not be a hard criteria in a zeolitic medium as 
there are well defined pores. Thus in zeolitic medium as long as the solutes 
are smaller than the pore size, the non monotonic size dependence can be 
observed when there is strong interaction between the zeolitic wall and the 
solute particles. We also note that for system 1, 2 and 3 the primary contri- 
bution to the diffusion arises from the binary collisions. In a porous zeolite 
medium we expect no density relaxation of the medium. This should lead to 
a complete decoupling between the solute dynamics and the medium. Thus 
ideally the collisions between the solute and the zeolite structure should de- 
termine the friction/diffusion. Hence, system 1,2 and 3 in some sense mimics 
the diffusion through a porous medium. 

The study further reveals that there is a strong one to one correlation 
between the non monotonic size dependence of diffusion and non monotonic 
size dependence of the peak value of the solute-solvent RDF. The RDF for 
system 1, where this anomalous behaviour for diffusion is observed, shows 
that for larger solutes the 1st peak is diffused but for smaller solutes the 1st 
peak is well defined. The solute which has the smallest value for the 1st 
peak of RDF shows the highest diffusion value. For other systems which do 
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not show any anomaly the smaller particles always have more diffused 1st 
peak (system 2 and 3) or the RDF has a weak dependence on the solute size 
(system 4). Thus we find that the effect of levitation is manifested in the 
RDF. 

Although we use fully self consistent MCT but for smaller solutes we find 
that the self consistent calculation for the self dynamic structure factor pre- 
dicts a long time tail which leads to larger friction value. The approximate 
self dynamic structure factor, using simulated mean square displacement and 
non-gaussian parameter does predict a fairly correct value for the diffusion. 
However for system 1 for the smallest two sizes there should be contribution 
from non hydrodynamic fast modes. Further, a more accurate memory func- 
tion for the self dynamic structure factor is required when the solute sizes 
are smaller. 

In conclusion, for the first time the anomalous size dependence of diffusion 
has been predicted using analytical MCT. The study reveals that for solute- 
solvent systems, this anomalous behaviour can be observed only when the 
solute-solvent interactions are strong and also the solute can explore the 
inter-solvent cage. The fact that the effect of a dynamic phenomena like 
levitation has its manifestation in the static RDF is quite encouraging. In 
future we would like to use analytical MCT to study diffusion in zeolitic 
medium. 
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